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1. Introduction 

Large progress in synthesis of heavy nuclei done in 
the last decades raises new demands for a better and 
better theoretical description of their decay modes. 
Spontaneous fission and o-radioactivity are the most 
important processes of disintegration of heavy nuclei. 

Alpha decay occurs most often in actinides region, 
but is also observed in isotopes with Z > 52. The 
first theoretical interpretation of alpha-decay process 
was given independently by Gamow [1], Gurney and 
Gondon [2] in 1928 year. Emission of alpha particle 
is treated as quantum-mechanical tunnelling through 
the nuclear Goulomb barrier, where the probability 
of emission is calculated using one-dimensional WKB 
approximation. Recently it was proven, that this 
approximation can be successfully used to evaluate the 
probability of tunnelling by alpha particle and cluster 
as well (see Ref. [3]). It was also shown, that within this 
simple model (containing only 1 adjustable parameter 
for even-even nuclei) one can reproduce alpha decay 
half-lives of heavy emitters with higher accuracy in 
comparison with modern (containing 5 parameters) 
version [4] of Viola-Seaborg formula [5]. 

The observation in 1938 of the neutron induced 
nuclear fission by Hahn and Strassmann came rather 
unexpected [6]. This new phenomenon was explained 
within a few weeks by Meitner and Frisch’s who 
established the most important features of low-energy 
fission: the energy released in this process is equal 
almost 200 MeV and results from the Goulomb 
repulsion of the fission fragments and the number 
of neutron emitted per fission event larger than one, 
what has opened a possibility for a chain reaction [7]. 
One and a half year after Hahn and Strassmann’s 
discovery, Flerov and Petrzhak have first observed 
the spontaneous fission of Uranium [8]. As we have 
written above, Gamow had explained the a-decay and 
its sometimes rather long half-lives by a quantum 
tunnelling process of a preformed a particle through 
a Goulomb barrier. So, according to the concepts 
of Meitner and Frisch, one could expect spontaneous 
fission of uranium also from the ground-state, but 
with a considerably longer half-life than for the a- 
decay, because of the larger reduced mass for almost 
symmetric fission. 

The first quantitative estimates of the sponta¬ 
neous fission probability becomes possible, when the 
microscopic-macroscopic model of the potential energy 


of deformed nuclei was developed and the Inglis crank¬ 
ing model was implemented to evaluate the inertia cor¬ 
responding to the fission mode. Of course in time these 
models give more and more precise reproduction of the 
spontaneous fission half-lives systematics and the fis¬ 
sion mass distributions. A short critical description of 
the spontaneous fission theories will be presented in 
our paper in order to better understand physical back¬ 
ground of a simple model for the probability of this 
decay mode, which we have developed following an old 
idea of W. J. Swiatecki. Namely, in 1955 year he pro¬ 
posed a formula which was able to describe the global 
systematics of the spontaneous fission half-lives [9]. 
His simple phenomenological formula, based on cor¬ 
relations between logarithms of observed spontaneous 
fission half-lives and ground state microscopic correc¬ 
tions, reproduced well the experimental data known 
at that time. We are going to show in the follow¬ 
ing that his idea, combined with the modern version 
of the liquid drop model (LSD) [10], allows to obtain 
satisfactory agreement with the measured up to now 
spontaneous fission half-lives. 

The paper is organized as follows: the main 
assumptions of a-decay model and results obtained for 
emitters with 52 < Z < 110 are presented in Sec. 2. 
The WKB theory of the fission barrier penetration as 
well as the the semi-empirical Swiatecki’s formula for 
spontaneous fission half-lives will be described in Sec. 
3, where results for the isotopes with 90 < Z < 114 are 
analyzed. Sec. 4. contains Summary. 


2. Alpha decay 


2.1. The model 


The quantum tunnelling theory of alpha emission 
assumes, that the decay constant A is proportional 
to the barrier penetration probability P, frequency of 
assaults on the nuclear Goulomb barrier per time-unit 
n and particle preformation factor Sa • In the presented 
model the expression for decay constant is simplified 
as the effect of the preformation is included into the 
probability P (see discussion in Ref. [3]): 

X = nP. (1) 


The probability P of tunnelling through the barrier 
is calculated using one-dimensional WKB approxima¬ 
tion: 


P = exp 


1 VMV{r) - E^)dr 


( 2 ) 
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Figure 1. Differences of the theoretical and experimental [11] Figure 2. The same as in Fig. 1, but for isotopes with 

alpha decay half-lives on logarithmic scale, calculated for nuclei 81 < Z < 90 (top) and 91 < Z < 101 (bottom), 

with 52 < Z < 70 (top) and 71 < Z < 80 (bottom). 
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differs from that reported in Ref. [3], as now more 
alpha emitters were taken and the cluster decays were 
not included in our analysis. Fitting procedures of 
hindmnee faetor were performed for odd systems with 
fixed ro value {h = 0.25 for odd-A and doubled for 
odd-odd emitters). 

2.2. Results 

The logarithms of the ratios of the half-lives 
calculated using above formalism, to the measured ones 
(^ 1 / 2 ) examined nuclei are shown in Figs. 1-3 

as a function of the neutron number N. 

Large deviations from the data are observed for 
emitters with neutron number N = 127 in gsBi 
- 92 U elements (Fig. 2), where underestimations of 
the logic(T 1 / 2 ) reach about two orders of magnitude. 
Calculated half-lives of some lighter emitters (52 Te - 
74W) are overestimated, but these discrepancies do not 
exceed one order of magnitude. The root-mean-square 
deviations of our estimations made for the a-decay 
half-lives for all considered nuclei are summarized in 
Table 1. 


Figure 3. The same as in Fig. 1, but for isotopes with 
102 < Z < 110. 


Table 1. Root-mean-square deviations of calcu¬ 

lated using Eq. (6) with the radius constant ro=1.23 fm. 


where /i is a reduced mass, R is the spherical square 
well radius 


R = ro{Al/^ + 


(3) 


and b is the exit point from the Coulomb barrier: 

6 =^= 1 ^. ( 4 ) 

and are the mass and proton numbers 

of an alpha particle (or cluster) and a daughter nucleus 
respectively. is the kinetic energy of emitted alpha 
(cluster) particle [11]. The number of assaults per 
time-unit n is evaluated from the quantum-mechanical 
ground-state frequency in the spherical square well: 


In this formalism the a-decay half-life can be expressed 
as follows: 


Ti/2 — 


ln 2 . 

T 


-lO"^ 


( 6 ) 


where constant h (so-called hindranee faetor) was 
additionally introduced for odd nuclei. The least- 
square fit of the radius constant ro was performed to 
known experimental half-lives of even-even a-emitters 
(127 cases) [ 11 ]. Only the most probable alpha decays 
mode at each isotope was chosen in our analysis. 
Obtained in this way value of the nuclear well radius 
constant is equal to ro = 1.23 fm, and slightly 


TTZ — TTat 

n 

h 

r.m.s. 

e-e 

127 

0 

0.39 

e-o 

97 

0.25 

0.66 

o-e 

82 

0.25 

0.54 

0-0 

54 

0.5 

0.79 


3. Spontaneous fission within the WKB theory 

A continuous interest in the theoretical description 
of the fission dynamics is observed since discovery of 
this phenomenon in 1938. In present Section we are 
going to concentrate on the spontaneous fission and 
its description within the WKB theory only. The 
results presented here will be used to understand, why 
a simple one parameter (for even-even nuclei) model 
a la Swiqteeki [9], described in Section 4, is able 
to reproduce the spontaneous fission half-lives of all 
known nuclei with higher accuracy than obtained using 
more advanced theories. 

Quantum mechanically it is possible for a nucleus 
in its ground state to tunnel the fission barrier. 
The probability of the barrier penetration depends 
not only on the height and width of the barrier 
but also on the magnitude of the collective inertia 
associated with the fission mode. The calculations of 
the spontaneous fission half-lives Tgf require a careful 
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evaluation of the collective potential energy surface and 
the collective inertia tensor. Different models, like the 
macroscopic-microscopic model [12, 13] or the Hartree- 
Fock-Bogolubov self-consistent theories (see e.g. [14]), 
can be used to generate the potential energy surfaces. 
The collective mass parameters one evaluates usually 
in the cranking approximation [15, 13] or within 
the Generator Coordinate Method (GCM) with the 
generalized Gaussian Overlap Approximation (GOA) 
[ 16 ]. 



Figure 4. A contour plot of the energy surface for is 

shown in the two-dimensional quadrupole (e) and hexadecapole 
( 64 ) deformation parameter space. The dashed line is the path 
of steepest descent from the saddle points to the ground state 
and to the scission line, the thick line is the minimal-action path. 
After Ref. [17]. 


rpsf _ 

^ 1/2 — 


( 7 ) 


The spontaneous fission half-life is given by: 

In 2 
nP 

Here n is the number of assaults per time unit, 
associated with the frequency of zero-point vibration 
of the nucleus in the fission mode direction. The fission 
barrier penetration probability P might be evaluated 
within the one-dimensional WKB approximation, what 
leads to the following expression: 

^ (8) 


P = 


1 + exp{25'(L)} 

The action-integral S{L), calculated along a fission 
path L{s) in the multi-dimensional space of collective 
coordinates is given by: 


S{L) = 


fl 




BM[V{s)-E^ 


1/2 


ds , 


( 9 ) 


where the integration limits Si and 52 , correspond to 
the classical turning points. is the energy of a 

nucleus in its ground state and V{s) is the collective 
potential. Bss{s) is an effective inertia tensor in the 
multi-dimensional space of collective coordinates {qi}: 

nn'i 

^ ds ds • ^ ^ 

k,l 


Bss{s) = Bki 




Figure 5. Collective potential (top) and inertia (bottom) 
along the static (thick solid line) and the dynamic (thin solid 
line) path to fission as a function of the relative distance between 
the fragment mass centres s = Ri 2 /Ro- The reduced (p,), 
irrotational fiow and phenomenological inertias [19] are shown 
for comparison. After Ref. [20]. 


Determination of the fission path in the deforma¬ 
tion space plays crucial role during evaluation of the 
fission probability [13]. The action integral (9) should 
be calculated along such a path (so-called dynamie)^ 
which minimize its total value, thereby maximize fis¬ 
sion probability (8) [18]. In contrast to the statie path^ 
the dynamic one (dependent on the inertia tensor) does 
not have to lead through the bottom of the fission val¬ 
ley, as one can see in Fig. 4. The comparison of the fis¬ 
sion barrier and effective inertia, corresponding to the 
dynamic and static fission paths for ^^^Fm is shown in 
Fig. 5. The collective potential V{s), corresponding to 
the dynamic path, is larger than in case of the static 
one, but the dynamic inertia is smaller. It might be 
roughly approximated by the so-called phenomenolog- 
ieal inertia fitted to the observed spontaneous fission 
half-lives [19]: 


M 


phen 

R 12 R 12 


17 

M ( 1 + ^^5 



( 11 ) 


Here r is the reduced mass, corresponding to the 
relative distance between fragments (i?i2=3/4 for a 
sphere) and a = 2.452 is a numerical constant. 
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N 

Figure 6. Spontaneous fission half-lives of even-even 
transactinides and their shape isomers (II) obtained with the 
WKB approximation (8) along the least action (dynamical 
trajectory). After Ref. [18]. 

determined from a fit to the exact irrotational flow 
inertia [21] for k = 1. It was found in Ref. [19] that the 
phenomenological inertia obtained with the constant 
k = 11.5 reproduces well the systematics of the known 
spontaneous fission half-lives. 

Spontaneous fission half-lives of even-even trans¬ 
actinides evaluated in Ref. [18] are presented in Fig. 6. 
The estimates correspond to the least-action trajecto¬ 
ries (dynamic path to fission). The WKB approxima¬ 
tion (8) and the cranking inertia tensor were used. The 
two-dimensional potential energy surfaces for each iso¬ 
tope were obtained within the macroscopic-microscopic 
model using the Nilsson potential and the droplet en¬ 
ergy for the macroscopic part. The effects of the axial 
asymmetry as well as the reflection asymmetry were 
taken into account. 

Usually calculations of spontaneous fission half- 
lives are performed with an assumption, that the 
collective potential is equal to the Hartree-Fock- 
Bogolubov (HFB) binding energy or its macroscopic- 
microscopic approximation. The zero-point energy 
corrections are often neglected in papers dealing with 
similar problems. We shell discuss this point later. 
Moreover, one assumes, that the ground state is located 
at AE ^ 0.5 MeV above the minimum of the HFB 
potential. As a justification of such a choice of 
AE one refers to experimental values of the average 


quadrupole phonon energy, which is about 1 MeV for 
actinides. Several authors even treat AE as a free 
parameter. It should be stressed, that the HFB theory, 
as each variational method, gives the upper limit of 
the ground-state energy, so the zero-point correction 
energy at the equilibrium point is approximately equal 
to the difference of the ground-state and the potential 
energy at the minimum (see Ref. [23]). It can be 
easily shown, that the zero-point energy of harmonic 
oscillations is equal to the half of the corresponding 
phonon energy, i.e. {^huj) [24], what is consistent 
with result, obtained in Ref. [25], where the coupled 
quadrupole and octupole vibrations in the region of 
Ra-Th nuclei were analyzed. This conclusion is also 
valid even for strongly anharmonic motion such as 
pairing vibrations [26]. 

The number (n) of the collective degrees of 
freedom results much on the value of the zero-point 
energy correction. Namely, the eo is proportional to 
U'huj. Thus, when the HFB equations are solved in the 
whole n-dimensional space of collective coordinates, all 
zero-point corrections should be substracted from the 
energy to obtain the collective potential: 

V{s) = Eufb{s) -5o®®(5) , (12) 

where is the zero-point energy related to the fission 
mode [16]. 

In the case of the one dimensional fission barrier 
penetration given by the Eq. (9), the oscillations 
perpendicular to the fission path are not included 
during evaluation of tunnelling probability (16). 
In this direction the zero-point energy corrections 
will be approximately cancelled by the 1/2 of the 
corresponding phonon energy. 

In Fig. 7 the potential and the binding energy 
along the fission path are plotted. As one can observe, 
the zero-point corrections may affect on the fission 
barrier height. 



Figure 7. Collective potential V{s), binding energy EhpbO) 
and ground-state energy Egs of a nucleus along the fission path 
s. After Ref. [20]. 

The determination of the least action path in 
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deformation space is not a trivial task (see Ref. [13] for 
further discussion). The approach, which is worth to be 
mentioned in this context is dynamical programming 
method [17, 18]. Nevertheless, applying this method 
one should choose optimal distance between the mesh 
points to avoid large numerical errors. That is because 
the accuracy of calculations of the partial derivatives 
^, which result much on the effective mass parameter 
(10) as well as the zero-point energy correction (see 
e.g. Ref. [16]) along the fission path, strongly depends 
on the grid size. The von Ritz implementation (see 
appendix of Ref. [27]) might be useful as an alternative 
for practical calculations of the spontaneous fission 
probability. The main advantage of this method is that 
the influence of the zero-point energy correction along 
the dynamical paths might easy be taken into account. 

Typically the mass yield of the spontaneous fission 
fragments of the actinide nuclei is spread between 
70 and 180. Usually, except the bimodal fission of 
258Fm and neighbouring nuclei, the asymmetric fission 
is more probable than the symmetric one and the most 
populated mass of the heavier fragment is around 140. 
A very mass asymmetric spontaneous fission process, 
when the mass of the lighter fragment is around 20, is 
already considered as a cluster emission. This type of 
radioactivity was predicted in 1980 by Sandulescu and 
co-workers [28] and four years later it was discovered 
by Rose and Jones [29], who observed the spontaneous 
emission of from ^^^Ra. The cluster emission is a 
very rare process. Its relative branching ratio to the 
a-decay is of the order 10“^^ to 10“^^. Nevertheless, 
in last two decades one has observed clusters from 
to ^^Si emitted by actinide nuclei from ^^^Fr to ^^^Cm. 
In all cases the residual nucleus was always close to the 
double magic ^^^Pb. 

230Th 


Figure 8. The HFB potential energy of 230 function of 
quadrupole and octupole moments. The the thick solid line show 
the path to cluster emission (CE) while the dashed one is drawn 
along the valley to asymmetric fission (AF). After Ref. [30]. 

In general, there are two alternative theoretical 


descriptions of the spontaneous cluster decay. One 
can use a fission-like mechanism to reproduce the 
main features of the process [28, 31, 32, 33, 30] or 
one assumes, in close analogy to the a-decay, that a 
cluster, performed by a nonadiabatic mechanism inside 
nucleus, penetrates the potential barrier created by 
the Coulomb and nuclear interaction with a daughter 
nucleus [34, 35, 36, 37]. The potential energy landscape 
of 230 Fp obtained by the HFB calculation with the 
Gogny DIS force [30] is presented in Fig. 8 as a function 
of quadrupole (Q 2 ) and octupole (Qs) moments. Two 
valleys leading to the cluster emission (CE, solid line) 
and to the asymmetric fission (AF, dashed line) are 
visible in the plot. 

The above theoretical models, which describe the 
cluster-radioactivity, are rather complex and contain 
adjustable parameters. Last year we have shown 
in Ref. [3], that within a simple Gamow model 
with only one adjustable parameter (radius constant) 
common for the a-decay and the cluster emission one 
can describe with good accuracy the experimental 
systematics of half-lives for the cluster radioactivity 
of even-even nuclei. An additional hindrance constant 
was introduced in Ref. [3] to describe the cluster 
emission probability from odd-even, even-odd and odd- 
odd isotopes. 

4. Simple phenomenological formula for the 
spontaneous fission half-lives 

Encouraged by the good result for obtained in 
the Gamow theory, we are going in the following 
to describe in another simple model the spontaneous 
fission half-lives of transactinide nuclei. Namely, we 
shall adopt the Swiatecki idea from 1955 [9], who found 
a simple relation between the spontaneous fission half- 
lives and the experimental mass deviations from their 
liquid drop estimates. The crucial ingredient of the 
model is the liquid drop formula which is described in 
the next subsection. 

4 . 1 . The liquid drop model 

In our analysis we decided to use the modern 
version of the macroscopic-microscopic model [10]. 
The macroscopic part of this Lublin-Strasbourg Drop 
(LSD) mass formula (in MeV units) is as follows: 

MLSD(^,A^,clef) = 

7.289034 • Z -1- 8.071431 • N - 0.00001433 • 

- 15.4920(1 - l.SmiP)A 

+ 16.9707(1 - 2.2938/2)^2/3 Hsurf (def) (13) 

+ 3.8602(1 + 2.3764/2)^1/3 Hcur(def) 

+ 0.70978 Z2/A1/3 Hcoui(def) - 0.9181 Z‘^/A 

- 10exp(-4.2|/|) Hcong(def) + Eo-e • 
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where the odd-even energy term £^e-o is given in 
Ref. [38]. In Eq. (13) A = Z + N denotes the 
mass number, / = {N — Z)IA reduced isospin and 
^surf, ^cur, ^Coui ^ud B^onz relative to the 
sphere: surface, curvature. Coulomb and congruence 
(see Ref. [38]) energies. The parameters in the first 
and the last row in Eq. (13) are taken from Ref. [39], 
while the rest 8 parameter were fitted to the data. 

4 . 2 . Simple model a la Swiqteeki 

Eollowing the idea presented in Ref. [9], we are going 
to find an approximative functional dependence of the 
logarithms of spontaneous fission half-lives, corrected 
by mass shifts k5M^ on proton Z and neutron N 
numbers: 

f{Z, N) = \og^^[Tl%/y] + k5M{Z, N) . (14) 

The 5M is a ground-state microscopic energy, defined 
as: 

N) = M,^p{Z, N) - Mlsd(^, V 0), (15) 

where N) is a measured mass of isotope, taken 

from [40] and Mlsd(^ 5 ^, 0) were calculated using 
formula (13). 



Figure 9. Logarithms of the observed spontaneous fission half- 
lives corrected with masses ’’shifts” kSM as a function of proton 
number. 

Eitting procedures included 35 fissioning even- 
even isotopes {Z < 102) with measured masses and 
half-lives. The smooth dependence for even-even 
isotopes was achieved for factor k = 7.7MeV“^, as it is 
shown in Eig. 9. The curves, fitted for odd-^d and odd- 
odd nuclei are shifted by a constant- hindranee faetor. 
It is equal to h = 2.5 for odd-A isotopes and doubled 
for odd-odd systems. Eormula for spontaneous fission 


half-lives in years is given by: 

logio (^) = -4-1 ■ min(Z, 103) + 380.2 

{ 0 for even — even , (16) 

2.5 for odd-A , 

5 for odd — odd . 



Z^/A 



/A Z^/A 

Figure 10. Spontaneous fission half-lives of even-even (a) 
and odd (b) nuclei, calculated using formula (16) (triangles) in 
comparison to the experimental values [11] (dots). 

Spontaneous fission half-lives, calculated using 
formula Eq. (16) are presented and compared with the 
data in Eig. 10. Isotopes which for masses are still not 
measured are marked with open symbols. A very good 
accuracy was achieved in actinides region, where the 
root-mean-square deviation of logT^^ equal to 1.18 
and it grows to 1.57 when odd-A and odd-odd nuclei 
are taken into account. The r.m.s. deviation reaches 
1.96, when one includes to the analysis the super-heavy 
isotopes, for which the liquid drop barrier vanishes. 
The quality of reproduction of known life-times by this 
simple model is even better than that presented in 
Eig. 6 obtained in the macroscopic-microscopic model 
with the cranking inertia [18]. Also more advanced 
modern calculations based on the HEB theory with 











On spontaneous fission and a-decay half-lives of atomie nuelei 


9 


the Skyrme (see e.g. [41, 42] or Gogny force [43] 
did not bring better estimates for large sample of 
isotopes than the simple formula presented above. It 
raises a question: why the simple Swiatecki model 
works? To understand it, one has to recall Swiatecki’s 
topographical theorem and the results obtained using 
dynamic (least action) trajectories to fission described 
in Section 3. 

4 . 3 . The topogmphieal theorem 



Figure 11. Experimental fission barrier heights compared with 
their estimates evaluated according to the topographical theorem 
of Swiatecki and the LSD macroscopic energy. After Ref. [46]. 

It was shown in Refs. [10, 27, 44] that the 
LSD model (13), which parameters were fitted to 
the experimental ground-state masses only, is able 
to reproduce well the fission barrier heights of light, 
medium and heavy nuclei, when the microscopic part of 
the ground-state binding energy is taken into account. 
According to the topographical theorem, proposed by 
Myers and Swiatecki [38], the mass of a nucleus in a 
saddle-point is mostly determined by the macroscopic 
part of its binding energy. The influence of the 
shell effects on the saddle-point energy is rather weak 
and in the first approximation might be neglected. 
So, one can evaluate the fission barrier height as a 
difference between the macroscopic (here LSD) and the 
experimental mass: 

Vb{Z, N) = Mlsd(^, N, saddle) - (17) 

where Mlsd (^5 A", saddle) is a macroscopic part of LSD 
mass formula (13) taken in a saddle point. 

The comparison of the experimental fission barrier 
height of actinide nuclei with their estimates, evaluated 
according to the topographical theorem of Swiatecki 
using the LSD macroscopic energy, is presented in 
Fig. 11. The r.m.s. deviation of the estimates from 
the data for all 18 considered fission barriers heights 
is 0.31 MeV only, what proves, that the topographical 


theorem really works [46]. All estimates do not exceed 
0.67 MeV and lie below the experimental values (except 
^^^Cf), what give a place for not washed out shell 
effects around the saddle point energy. 

4.4- Justifieation of the simple formula for Tgf 

Let us consider one dimensional fission barrier along 
the least-action (dynamic) trajectory described in 
Sec. 3. Assuming, that the fission path is parametrised 
by the collective coordinate s, one can evaluate 
the potential V{s) and the mass parameter Bss{s) 
corresponding to this path. The fluctuations of the 
inertia Bss along the least-action trajectory becomes 
smaller and smaller, when one increase the number 
of collective coordinates. Especially the collective 
pairing degrees of freedom significantly wash out the 
fluctuations of the inertia function [47]. 



Figure 12. Schematic plot of a fission barrier in form of two 
inverted parabolas smoothly joined at the scission point. 


A simple transformation: 



"^sadd 


(18) 


bring us to the new coordinate x{s), in which the 
inertia becomes constant Bxx = m. The lower 
integration limit in Eq. (18) is chosen at point x = 0, 
corresponding to the saddle point Sgadd- The collective 
potential in the new coordinate, schematically plotted 
in Eig. 12, can be approximated by two inverted 
parabolas of different stiffnesses Ci and Cr having the 
maximum at the saddle point: 


V{x) = 


ksadd 2 ^ 


for 

for 


X < 0 , 
X > 0 , 


(19) 


The stiffnesses of the V potential are chosen in such a 
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way that the action-integral S (9) becomes equal: 

K - - 

I^J^[V{x) - Eo] dx=j - Eo] dx (20) 

-x[ -XI 

where pairs and {—xi^Xr) are classical left 

and right turning points for the true and approximative 
potential respectively. The last integral in Eq. (20) can 
be rewritten as 



where Vb = V^add ~ Eq is the fission barrier height. 
After a small algebra the action integral becomes 


5' = 




TT UJi + Ur 
fl 2uJiUJr ’ 


( 21 ) 


where uji = ^/ClJm and ujr = are frequencies 

of the left and right (inverted) oscillator respectively. 
Introducing the average oscillator frequency 


2ujiu)j> 
LUl + iUr 


( 22 ) 


one can bring the action integral to the following form: 


S = 


2Vb 

hu) 


(23) 


In our approximation the action integral is propor¬ 
tional to the fission barrier height measured in the en¬ 
ergy quanta of the harmonic oscillator, which approx¬ 
imates the fission barrier form. For the action-integral 
S > 1 the logarithm of the spontaneous fission half- 
lives (7) can written as 

logiTf/Jy) = 2S- log(n) - 0.3665 (24) 

where the constant 0.3665=log[ln(2)] and n is the 
number of assaults of nucleus on the fission barrier 
per year {y). Having in mind, that according to 
the topographical theorem the fission barrier height is 
Vb = ~ ^exp making use of the relation 

(23), one can rewrite the last equation as follows: 

log{TlL/y) + —?^ = —^-log(n)-0.3665 ,(25) 

' TkjO Tluj 

where was defined in Eq. (15) and = 

MllV - is the liquid drop fission barrier 

height. The right hand side of Eq. (25) is a very 

smooth function of nucleon numbers as it is defined 

only by global properties of nucleus. Note, that 
the derived equation has the same structure as the 
phenomenological formula of Swiatecki (14), what 
proves his Ansatz. 

The liquid drop barrier height of actinides 
decreases almost linearly in function of Z from 4.3 MeV 
for Z = 90 to 0 for Z > 103 [44]. The fission barrier of 
finite height appears in the super-heavy nuclei mostly 



Figure 13. Logarithms of the observed spontaneous fission half- 
lives corrected with 7.7SM as a function of liquid drop barrier 
height. 


due to the shell effects in the ground state. The smooth 
dependence of logarithms of spontaneous fission half- 
lives, corrected by the ground state shell-plus-pairing 
effects, on the LSD fission barrier heights, shown in 
Fig. 13 for even-even (e-e), odd A (o-A) and odd-odd 
(o-o) nuclei. The data for e-e isotopes lie very close 
to the straight line, what validates Eq. (25). The data 
for o-A and o-o are above this line, what is due to the 
specialization energy, which increase the fission barrier 
heights. As it was mentioned before, stability of super¬ 
heavy nuclei is determined by shell effects. The liquid 
drop barrier of these isotopes vanishes, what is visible 
in Fig. 13 for nuclei with Z > 103. 

5. Conclusions 

The following conclusions can be drawn from our in¬ 
vestigation: 

• Simple, consistent model was applied to reproduce 
alpha decay half-lives in 360 nuclei with atomic 
number 52 < Z < 110. 

• Model reproduces o—decay half-lives with quite 
good accuracy; the root-mean-square deviation of 
logio(T^^ 2 ) for even-even isotopes is equal to 0.39. 

• Large underestimations in half-lives of N = 127 
isotones arises from strong shell effects, which are 
not considered in this simple model. 

• Semi-empirical formula for the spontaneous fission 
half-lives, depending on proton number and the 
ground state microscopic corrections, reproduces 
data for even-even super-heavy nuclei with rea¬ 
sonable accuracy. 

• Quality of spontaneous fission half-lives evalua¬ 
tion breaks down for nuclei with not measured yet 
masses. 
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• The logarithms of the spontaneous fission half- 
lives, corrected by the ground state shell-plus- 
pairing effects, are roughly proportional to the 
macroscopic barrier heights in nuclei up to Z=102. 
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